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Brief 


An  extension  of  the  Generalized  Standard  Addition  Method  is 
made  to  include  usage  of  more  sensors  then  analytes  and  signal 
drift  detection  and  correction. 


Abstract 


The  detection  and  subsequent  correction  of  drift  in  an  analy¬ 
tical  signal  occurring  during  the  course  of  an  analysis  is  a  common 
problem  in  analytical  chemistry.  In  this  paper,  a  solution  to  this 
problem  is  presented  in  the  form  of  an  extension  to  the  theory  of  the 
Generalized  Standard  Addition  Method  (GSAM) ,  a  method  currently  used 
to  simultaneously  correct  for  matrix  effects  and  interferences  in 
multicomponent  analysis.  Additionally,  the  GSAM  has  been  extended 
to  detect  and  correct  for  normally  unwanted  interference  effects  by  using 
more  sensors  (wavelengths,  electrodes,  etc.)  than  there  are  analytes. 
Standard  additions  of  a  pseudo-analyte  containing  all  potential  interfering 
components  In  a  sample  can  be  used  to  correct  for  all  potential  inter¬ 
ferences  at  once. 
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INTRODUCTION 

Until  fairly  recently,  analytical  chemists  have  used  laboratory 
computers  primarily  for  data  acquisition,  storage  and  retrieval.  Presently, 
a  growing  awareness  and  application  of  such  topics  as  experimental  design, 
optimal  control  and  multivariate  data  analysis  promises  to  expand  the  role 
of  computers  in  chemical  analysis.  With  greater  utilization  of  the  com¬ 
putational  and  logic-based  decision-making  abilities  of  computers,  a  new 
generation  of  what  might  be  called  intelligent  analytical  instrumentation 
will  evolve  with  the  ability  to  detect  and  correct  for  problems  that  may 
render  an  analysis  invalid. 

One  well-known  problem  is  drift  in  an  analytical  signal  occurring 
over  a  period  of  time  in  which  a  sample  is  analyzed.  This  is  called  short 
term  drift  as  opposed  to  long  term  drift  which  occurs  from  sample  to  sample 
or  day  to  day.  When  the  simple  Standard  Addition  Method  or  the  Generalized 
Standard  Addition  Method  (1)  is  used,  long  term  drift  is  not  a  problem  as 
calibration  and  analysis  is  performed  simultaneously  for  each  sample.  Short 
term  drift  may  be  a  consequence  of  various  factors  —  temperature,  electronic 
drift  in  detectors,  changes  in  optical  alignment  and  other  causes  depending 
on  the  analytical  instrument  of  interest.  Without  proper  drift  correction 
In  quantitative  chemical  analysis,  a  time-dependent  systematic  (non-random) 
error  is  introduced.  To  compensate  for  drift,  prior  knowledge  of  its 
presence  must  be  available.  This  is  a  severe  limitation.  Most  certainly, 
a  method  is  needed  that  can  detect  the  presence  of  short  term  drift,  charac¬ 
terize  the  drift,  and  simultaneously  correct  for  it.  Ideally,  this  method 
should  be  applicable  to  multicomponent  chemical  analysis. 


Recently,  our  laboratory  has  developed  an  experimental  design  and  data 
analysis  procedure  that  is  capable  of  fulfilling  these  requirements  (1,2). 
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It  is  a  generalization  based  on  the  method  of  standard  additions  and  can 
include  any  number  of  analytes.  It  Is  referred  to  as  the  Generalized 
Standard  Addition  Method  (GSAM)  and  was  originally  developed  to  correct 
for  matrix  effects  and  interferences  (i.e.,  chemical,  physical,  and 
spectral  interferences)  arising  from  chemical  species  in  the  sample.  It 
has  been  applied  to  spectrophotometry  (2) ,  inductively  coupled  plasma- 
atomic  emission  spectrometry  (ICP-AES)  (3),  ion-selective  electrode  poten- 
tiometry  (4),  and  anodic  stripping  voltammetry  (5).  In  this  work,  the 
GSAM  is  extended  to  detect  and  correct  for  interferences  that  do  not 
originate  from  chemical  species  present  in  the  sample  (i.e.,  instrumental 
drift) . 

When  the  response  from  an  analytical  sensor  is  drifting  for  one  reason 
or  another,  the  response  is  a  function  of  the  concentrations  of  analytes 
and  interfering  components,  and  the  cause  of  the  drift  (time,  temperature, 
etc.).  Before  a  drift  correction  can  be  made,  the  drift  must  be  detected 
in  much  the  same  way  as  the  analyte  or  an  interferent  is  detected.  The 
problem  is  one  of  detecting  the  drift,  a  nonrandom  signal,  in  the  presence 
of  the  random  noise  inherent  to  the  sensor.  To  test  the  potential  of  the 
GSAM  to  solve  the  drift  problem,  it  was  necessary  to  use  computer  simulated 
experiments  generating  carefully  controlled  data  to  be  analyzed.  This  is 
exactly  analogous  to  preparing  primary  standard  samples  in  order  to  test 
the  recovery  or  detection  limits  of  a  new  analytical  method.  The  method 
used  to  simulate  the  multicomponent  analysis  and  the  results  obtained  from 
the  GSAM  are  given  in  a  later  section  along  with  results  using  experimen¬ 
tally  obtained  data. 


In  addition  to  the  use  of  the  GSAM  for  detection  and  correction  of 


short  term  drift,  this  paper  will  also  extend  the  theory  of  the  GSAM  for 
the  case  where  the  number  of  sensors  (p)  is  greater  than  the  number  of 
analytes  (r). 

A  common  problem  in  chemical  analysis  is  the  presence  of 
potential  interfering  species  not  expected  in  the  sample.  This  is  an 
ideal  case  for  using  more  sensors  than  analytes.  For  example,  in  the 
routine  analysis  of  two  analytes  using  the  GSAII,  standard  additions  of 
two  standards  corresponding  to  the  two  analytes  should  be  made.  But 
if  potential  interfering  chemical  species  are  present,  experimentally 
recovered  concentrations  could  be  in  error.  By  using  more  sensors  than 
analytes,  this  paper  will  show  how  the  GSAM  can  be  used  to  indicate  the 
presence  of  potential  interferants  and  correct  the  analyte  concentrations 
accordingly. 
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Background  and  Theory 
The  Generalized  Standard  Addition  Method 

The  theory  of  the  GSAM  has  been  presented  in  full  detail  in 
previous  papers  (1,2).  Therefore,  only  a  basic  review  is  given  here. 
The  GSAM  requires  that  for  r  analytes,  the  responses  from  p  sensors 
(p  £  r)  be  recorded  before  and  after  n  standard  additions  are  made 
(n  -  r) .  The  model  is 


I  c  k  o 

I*1  m,s  s,£ 

3~1 


m-1, . « . ,n 
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(1) 


where  r  „  is  the  response  of  the  Jlth  sensor  after  the  mth  addition 
m,  JC 

of  analyte  (or  interfering  component)  s.  c  is  the  concentration  of 

m,s 

the  sth  component  and  k  0  is  the  linear  response  constant  for  the  £th 

S  y  Jo 

sensor  to  the  sth  component.  Representing  the  above  in  matrix  notation 
gives 

R  =  CK  (2) 


where  R  is  the  (n  +  1)  x  p  matrix  of  measured  responses  (initial 
responses  plus  the  responses  following  the  n  additions) ,  C  is  the 
(n  +  1)  x  r  concentration  matrix,  and  K  represents  the  r  x  p  matrix  of 
linear  response  constants  showing  the  contribution  of  each  of  the  r 
analytes  to  each  of  the  p  sensors  for  a  given  sample  matrix. 

Since  concentrations  are  not  always  additive,  volume  corrected 
responses  are  used.  Equation  (1)  can  be  rewritten  as 


lm,£ 


r  n 
r  mAs 

i  v 

s=l  m 


(3) 


where  n  is  the  total  number  of  moles  of  component  s,  and  V  is 
m,s  r  m 

the  volume  after  the  mth  addition.  This  can  be  expressed  as 
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a  ,  =  Vr  .  *  J  n  k  „ 
m  m.i,  L,  m.s  s.x. 
s=l 


(4) 


which  can  be  separated  into 


Sn.S. 


An 


(5) 


where  An  is  the  total  amount  of  component  s  added  and  n 

™>s  o,s 

initial  number  of  moles  of  component  s.  Now 


is  the 


1. 

1  a  =  V  r  =  y  n 
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(6) 


with  Vq  representing  the  initial  volume  of  the  sample  before  any 

additions  are  made  and  r  „  the  initial  response  of  the  Jlth  sensor. 

o,l 

Equation  (5)  can  now  be  written  as 


q  =  J  An  k  „  +  q  „ 
m,£.  ^  m.s  s,Jt  o,£ 


Considering  only  the  changes  in  responses  gives. 


\q_  ,  ■  I  An  k  „ 
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In  matrix  notation. 


AQ  -  ANK 


(7) 


(8) 


(9) 


where  AQ  is  nxp,  AN  is  n  x  r,  and  k  is  r  x  p.  Using  the  method  of  multiple 
linear  least  squares  to  solve  for  K  in  the  presence  of  an  overdetermined 
system  (n>r) ,  the  generalized  inverse  of  AN  (6)  is  needed  resulting  in, 

K  -  (ANTAN)_1  ANTAQ  (10) 

where  the  generalized  inverse  of  the  matrix  A?* is  defined  as  (ANTAN)_1AnT 


After  K  is  found,  the  vector  of  initial  analyte  quantities,  nQ  , 
is  recovered  by  solving  (when  r  =  p) 


n 

o 


(ID 


Therefore,  to  perform  an  analysis  using  the  GSAM,  three  conditions 
must  be  met: 

(1)  n  -  r,  (2)  the  standard  additions  must  span  the  r  dimensional 

T 

concentration  space  (AN  AN  must  have  rank  r) ,  and  (3)  each  analyte  must 

T 

effect  a  response  from  at  least  one  sensor  (AQ  A Q  must  have  rank  r). 

The  GSAM  is  a  very  simple  method  to  employ.  To  determine  one  row 
of  an  r  x  p  K  matrix  (the  sensitivities  of  the  p  sensors  to  the  sth 
component),  only  one  addition  containing  the  analyte  is  required  while 
recording  all  the  sensor's  responses.  Equation  10  then  simplifies  to 


_1_ 

An 


(12) 


This  process  can  be  repeated  n  =  r  times  making  standard  additions  of 
each  analyte  only  once  permitting  the  determination  of  K  and  thereby 
allowing  nQ  to  be  estimated.  Unfortunately,  utilizing  this  simple 
approach  the  analyst  has  obtained  just  sufficient  data  to  satisfy  the  above 
conditions.  With  n  >  r  a  better  estimation  of  the  elements  of  K  is  obtained 
due  to  the  use  of  an  overdetermined  system  of  equations.  Nevertheless, 
the  simplified  mode  of  the  GSAM  permits  a  very  quick  estimate  of  inter¬ 
ferences  present  and  initial  amounts. 

Drift  Correction 

Only  small  alterations  in  the  theory  of  the  GSAM  ara  needed  to  yield 


a  model  capable  of  compensating  for  drift.  For  first  order  drift 


equation  (3)  becomes 


r  n 


=  V  JUS.  k  +  t  k 
m,2.  ^  V  s,t  m  t,J 

s=l  m 


(13) 


where  t  is  the  time  that  has  elapsed  from  the  initial  measurement  up 


m 


to  the  mth  addition  and  k  .  is  the  linear  response  constant  of  time 

t  »*> 

on  the  Slth  sensor.  Equation  (4)  becomes 


L 

a  „  =  7  n  k  .  +  V  t  k  „ 

Tn,ji  ^  m,s  s,x.  m  m  t.x. 


(14) 


and  equation  (6)  is  modified  to  give 

r 
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o,Z 


I 

s=l 


n 

o ,  s 


k  .  +  V  t  k 
s,x.  o  o  t 


Z 


(15) 


where  t  is  the  initial  time.  Usually,  the  last  term  will  be  zero 
as  tQ  will  be  started  at  zero.  Equation  (14)  can  be  expanded  to  give 


r  r 

7  An  k  ,.  +  )  n 

,  m,s  s ,Z  ,  o,s 
s=l  s=l 


k  o  +  v  t  k 
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(16) 


Subtracting  equation  (15)  from  equation  (16)  results  in 


Aq  =  y  An 
m,£.  L.  m, 
s=l 


s  ^s.X. 
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Using  matrix  notation  equation  (9)  is  again  obtained  except  in  this 

case  AN  is  now  n  x  (r  +  1)  where  the  entries  in  the  (r  +  l)th  column  are 

(V  t  -  V  t  )  =  At  ,,  for  the  mth  addition  and  K  becomes  (r  +  1)  x  p 
mm  oo  m,r+l  K 

with  the  (r+l)th  row  estimating  the  effect  of  time  on  each  of  the  p 
sensors.  In  the  general  case,  the  above  equations  can  be  mouifieu 

to  include  combinations  of  higher  order  functions  of  time.  For  these 
instances,  the  model  (eq.(l))  becomes 
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m,  H 


r 
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3=1  m*S 


"s,  Z 


w 

I 

i-=l 


(18) 


where  w  represents  the  order  of  the  drift  process.  t  should  be  noted 

that  the  model  can  be  modified  to  include  other  parameters  that  could 

cause  drift,  i.e.  temperature  cycling. 

Equation  (10)  is  still  used  to  solve  for  K.  However,  since  there  is 

no  time  sensor,  the  respective  rows  of  K  containing  response  constants 

for  time  must  be  deleted  before  using  equation  (11)  to  solve  for  n  . 

o 

These  deletions  are  allowed  because  by  making  the  "time  additions"  (including  time 
in  the  model) ,  K  is  corrected  for  drift  and  there  is  no  need  of  the  ktv 
rows.  Nevertheless,  ktW  is  quite  important.  In  view  of  the  fact  that 

ktW  represents  drift,  inspection  of  ktW  will  reveal  to  the  analyst  whether 
or  not  drift  is  present  in  any  of  the  p  sensors.  Student  t-tests  can  be 
performed  to  detect  if  a  specific  ktW  ^  is  significantly  different  from 
zero  implying  the  presence  of  drift.  Thus,  the  GSAM  will  inform  the 
operator  when  drift  is  present  and  concurrently  correct  for  it.  When  used 
in  an  automated  system,  the  computer  can  intelligently  correct  for  drift 
and  flag  the  sensor  as  potentially  defective. 

With  the  inclusion  of  the  so-called  time  additions,  the  three  guide¬ 
lines  mentioned  earlier  become  (1)  nk  (r+  number  of  time  parameters). 

T 

(2)  AN  AN  must  have  rank  (r  +  number  of  time  parameters) ,  and  (3)  remains 
T 

the  same  (A Q  AQ  must  have  rank  r) 

In  a  recent  publication  (2),  different  variations  of  the  GSAM  were 
introduced.  Due  to  the  consumption  of  sample  inherent  with  some  instrumentation, 
(i.e.  IC3-AES),  and  for  reasons  to  be  considered  later  in  the  paper,  the 
incremental  difference  computation  partition  GSAM  (IDC-PGSAM)  was  the 


most  appropriate  variant  to  use. 


9 


The  IDC  dictates  that  the  entries  of  AN  and  Atyare  just  the  amounts  added 
or  time  elapsed  since  the  previous  addition.  Likewise,  the  AQ  entries 
for  IDG-PGSAM  are  the  incremental  differences  in  the  volume  corrected 
response  from  the  previous  responses.  The  PGSAM  consists  of  partitioning 
the  sample  into  r  aliquots  with  additions  of  only  one  standard  for  each 
aliquot. 

More  Sensors  than  Analytes  (p  >  r) 

To  use  more  sensors  than  analytes,  the  GSAM  requires  trivial  alter¬ 
ations  in  the  equations.  Equations  (1)  through  (10)  remain  the  same. 
Solving  for  n^,  however,  requires  a  new  procedure.  Since  p  >  r,  K  is  not 
a  square  matrix  and  consequently  K  ^  no  longer  exists  invalidating 

equation  (11).  Therefore  a  generalized  inverse  of  K  is  needed  which  will 
the  estimatation  of  nQ  as 


T  T 
n  =  q  K  [KK1] 
o  o 


-1 


(19) 


T  T  —1 

where  the  generalized  inverse  of  the  K  matrix  is  defined  as  K  (KK  )  . 

The  guidelines  on  the  number,  composition,  and  order  of  the  n  multiple 


standard  additions  become  (1)  remains  the  same  (n  >  r) ,  (2) 
(ANlAN  must  have  rank  r) ,  and  (3)  AQTAQ  must  have  rank  £  r. 


The  IDC-GSAM  was  the  variation  used  in  the  p  >  r  study  (all  standard 
additions  are  made  to  the  same  sample) . 


allow 


Detection  of  Potential  Interferences 

The  GSAM  will  correct  for  the  interferences  during  a  multicomponent 
analysis  only  if  the  interferences  are  anticipated  prior  to  analysis. 
Correction  for  totally  unaccounted  for  interferences  is  considerably 
more  complicated  depending  upon  how  the  interfering  species  affects  the 
sensors . 

If  the  goal  of  a  routine  analytical  procedure  is  the  analysis  of 
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r  analytes  and  U  additional  potential  interferents  can  be  identified, 
then  the  GSAM  can  be  modified  to  detect  and  correct  for  the  presence 
of  one  or  more  of  the  U  potential  interferents  during  the  analysis  of 
the  r  analytes. 

Major  modifications  are  not  necessary  in  order  to  incorporate  this 
new  improvement  into  the  GSAM.  Equations  (1)  through  (11)  essentially 
remain  the  same.  The  main  alterations  required  lie  in  the  experimental 
design.  For  instance,  if  two  analytes  are  involved  (r=  2),  then  two  sensors 
(p  =  2)  and  two  or  more  standard  additions  (n  £  2)  would  normally  be 
chosen.  For  detection  and  correction  of  unaccounted  potential  interferents 
present  in  the  sample,  additional  sensors  yielding  responses  to  the 
potential  interferents  are  obviously  required  along  with  at  least  one 
extra  standard  addition  consisting  of  a  pseudo-analyte.  This  pseudo-analyte 
should  consist  of  one  solution  containing  all  of  the  U  potential  interferents. 
The  responses  for  all  sensors  are  recorded  as  usual.  Averaging  the 
extra  responses  for  the  potential  interferents  to  form  one  pseudo-reponse 
and  therefore  one  pseudo-sensor,  and  making  only  one  standard  addition 
containing  all  the  U  anticipated  interferences,  the  dimension  of  AQ,  AN,  K  be¬ 
come  respectively,  n  x  (p+1) ,  n  x  (r+1)  and  (r+1)  x  (p+1).  Provided  additions 
of  all  potential  analytes  are  made,  the  estimated  n^  will  be  interference 
corrected  with  the  moles  of  pseudo-analyte  corresponding  to  the  average  moles 
of  all  potential  interferents.  Naturally,  when  none  of  the  U  interferents 
are  present,  the  GSAM  estimated  moles  for  the  pseudo-analyte  will  be  zero. 

If  one  or  more  potential  interferents  are  present,  the  pseudo-analyte  will 
be  non-zero  and  the  primary  analytes  will  be  correctly  determined.  This 
procedure  makes  no  provision  for  determining  exactly  which  of  the 
pseudo-analyte  inter ferent  is  actually  present  in  the  sample.  Only  the 
primary  analytes  are  of  interest. 


11 


EXPERIMENTAL 

All  calculations  were  done  using  "GSAM",  a  Fortran  IV  program 
available  from  Infometrix,  Inc.,  P.0.  Box  25808,  Seattle,  WA  98125.  All 
computations  were  performed  on  the  Department  of  Chemistry,  VAX  11/780 
computer. 

Time  Drift  Study 

All  "synthetic"  data  from  computer  simulated  analyses  (aq  as  expressed 

in  eq.  (9))  were  based  on  known  K  and  AN  matrices.  Three  analytes 

were  examined  with  three  standard  additions  of  1  unit  each  per  analyte. 

A  second  order  time  model  was  used  In  this  study.  Linear  time  additions 

consisted  of  1  unit  of  time  (units  are  arbitrary)  between  additions  of 

analyte  in  a  given  partition  with  2  units  of  time  elapsing  between  partitons. 

Quadratic  time  additions  were  performed  in  a  similar  manner.  Volumes 

were  kept  constant.  A  2%  quadratic  drift  was  incorporated  into  the 

synthetic  data  responses  (no  random  noise)  of  sensor  three  corre¬ 
sponding  to  analyte  three  and  a  2%  linear  drift  was  simulated  in  sensor 

one  corresponding  to  analyte  1.  Random  noise  was  introduced  to  the  simu¬ 
lated  data  responses  by  using  a  Monte  Carlo  method  (7).  A  normal  distribu¬ 
tion  was  assumed  with  a  mean  zero  and  a  standard  deviation  equal  to  a  1% 
relative  standard  deviation  (RSD)  of  the  noise-free  synthetic  data  responses. 
A  1%  RSD  is  the  estimated  precision  of  many  analytical  instruments.  Twenty 
Monte  Carlo  perturbations  were  performed  to  compute  twenty  random 
responses.  From  the  final  results  using  the  twenty  Monte  Carlo 

responses,  standard  deviations  were  estimated  for  n  and  K. 

o 

~~7 

Table  I  lists  the  synthetic  data  before  introduction  of  random  noise.  As 
mentioned  earlier,  the  IDC-PGSAM  experimental  design  was  used. 


The  ICP-AES  and  experimental  operating  conditions  used  have  been 
described  elsewhere  (3) .  Also  presented  in  this  reference  is  the  pro¬ 
cedure  of  additions  and  solution  concentrations.  A  brief  summary  of 
the  experimental  design  is  listed  in  Table  II.  Time  additions  consisted 
of  three  minutes  between  analyte  additions  and  five  minutes  between 
partitions.  A  10%  linear  drift  was  mathematically  added  to  the  A1  sensor. 

Standard  deviations  for  n  and  K  were  calculated  with  the  Monte  Carlo 

o 

method  using  twenty  perturbations  with  a  RSD  equal  to  1%,  the  estimated 
precision  of  our  ICP  instrument. 

Potential  Interferences  Study 

The  time  independent  data  (no  drift)  used  in  the  simulated  time  study 
(Table  I)  was  used.  Sensor  three  represented  the  pseudo-sensor.  Only  one 
standard  addition  composed  of  all  unforeseen  analytes  was  made.  Random 


responses  were  calculated  as  stated  above. 
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RESULTS  AND  DISCUSSION 
Time  Drift  Study 

Four  simulated  experiments  were  performed  in  order  to  test  the  ability 
of  the  GSAM  to  detect  and  correct  for  time  drift  while  avoiding  over  compen¬ 
sation  when  drift  is  not  present.  They  consisted  of:  no  time  additions  and 
no  drift;  time  additions  but  no  drift;  no  time  additions  and  with  drift  present; 
and  time  additions  with  drift.  Results  from  these  four  studies  are  presented 
in  rabies  III,  IV,  V  and  VI,  respectively.  Before  examining  the  results 
given  in  the  Tables,  some  clarification  of  the  experimental  procedure  is 
necessary.  From  the  experimental  section  it  is  noted  that  the  IDC-PGSAM 
was  used  and  the  corresponding  At  in  AN  are  incremental.  Namely,  the  time 
elapsed  since  the  previous  measurement  is  used  for  the  entries  in  the  Aty 
For  the  second  partition  the  first  linear  time  addition  (fourth  time  addition 
overall)  coinciding  to  the  first  standard  addition  of  the  same  partition  is 
five  units  of  time  and  not  the  expected  two  units.  This  is  true  since  the 
partitions  are  thought  of  as  new  samples  and  Atwin  this  case  is  the  time 
elapsed  since  the  initial  reading  was  recorded  (three  units  plus  two  units 
for  the  time  expired  between  partitions) .  The  change  in  the  increments  of 

time  between  partitions  is  necessary  to  preserve  the  non- singularity  of 
T  -1 

(AN  AN)  .  Analogous  inference  is  given  to  the  other  entries  of  the  At 

2 

corresponding  to  the  additional  partitions  and  to  At  .  Also,  while  monitor¬ 
ing  the  time  between  additions  it  is  assumed  that  during  the  brief  instant 
required  for  measurement  of  the  signal,  drift  is  negligible. 

From  Table  I  it  is  observed  that  no  units  of  time  have  been  ascribed 
to  Aft  Actually,  the  units  for  Atware  immaterial.  Units  of  time  only  alter  the 

magnitude  of  kt  and  kfc2  (last  two  rows  of  K,  effect  of  time  or  drift  on  each 
■7  ~r 

sensor)  and 
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not  the  's  which  remain  constant.  The  magnitude  of  the  entries  in 
At  or  At2  (decided  by  the  units  used)  will  then  dictate  the  magnitude 
of  drift  seen,  if  any,  in  the  sensors.  Therefore,  the  changing  or 

with  units  of  tine  are  of  no  consequence  since  k^  and  k^.2  are 

—7  -7  -? 

deleted  from  K  before  calculating  n  and  the  k  's  remain  time-corrected. 

o  r 

Finally,  from  Table  I  it  is  seen  that  an  additional  sensor  to 

account  for  time  additions  is  not  needed.  To  use  a  time  sensor  would  only 

allow  the  calculation  of  t  ,  but  t  is  already  known. 

The  results  in  Table  III  show  the  performance  of  a  normal  run  of  the 

GSAM,  that  is,  before  time  additions  are  made  and  drift  is  simulated.  The 

small  relative  errors  in  n  and  the  standard  deviations  of  n  indicate 

o  o 

that  the  GSAM  can  easily  tolerate  1%  random  noise  without  introducing 
large  errors  or  uncertainties  in  n  .  A  t-test  at  the  95%  confidence 
level,  shows  that  k2  ^  k3  ^  and  k1>3  (k^  column)  are  significantly 

different  from  zero  as  expected  according  to  Table  I. 

Table  IV  exhibits  the  results  when  there  is  uncertainty  about  the 

presence  of  drift.  That  is,  time  additions  are  made  without  prior 

knowledge  of  the  existence  of  drift.  Once  again,  using  t-test  (at  the  95% 

confidence  level)  will  show  that  k?  . ,  k.  ,  k  ,  k  ,  and  k  ?  are  not 

statistically  different  from  zero.  Also,  a  standard  t-test  at  the  95% 

confidence  level  indicated  that  there  is  no  statistical  difference  In  n 

o 

from  Tables  III  and  IV. 

Table  V  clearly  shows  that  the  errors  incurred  if  drift  is  present 

and  no  time  additions  are  made  to  compensate  for  the  drift  are  large 

indeed.  V.  ,  k  ,  and  k  are  now  significantly  different  from  zero. 

J  »  i- 

Although  the  rest  of  K  is  approximately  equal  to  the  previous  K,  K  is 
seen  to  be  substantially  dissimilar  to  K  ^  represented  in  Tables  II I  and  IV 
Accordingly,  since  K  ^  is  used  in  estimating  nQ,  the  errors  in  nQ  are 
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1 

i 

I 

i 

greatly  enhanced. 

Finally,  Table  VI  presents  the  results  when  time  additions  are 
made  and  drift  is  present.  k2  and  k^  3  are  once  again  not 

statistically  different  from  zero  and  the  GSAM  has  correctly  character¬ 
ized  the  magnitude  of  the  existing  drift.  Specifically,  k  and  k  2 

t , 1  t  ,  j 

are  approximately  equal  to  0.02.  K  ^  is  back  to  its  correct  form  (time- 
corrected)  and  n^  is  accurately  determined. 

From  the  results  presented  it  appears  the  GSAM  can  detect,  charac¬ 
terize,  and  correct  for  drift  in  the  presence  of  sample  related  interferences 
and  random  noise.  Other  studies  were  performed  to  determine  the  quantity 
of  random  noise  permissible  before  the  GSAM  would  fail.  With  the  K  in  Table  I, 
up  to  5%  random  noise  could  be  tolerated  before  relative  errors  in  n 

o 

approached  8%. 

However,  another  variant  of  the  GSAM  commonly  used  in  analysis  is 
the  IDC-GSAM  where  the  sample  is  no  longer  partitioned,  and  all  standard 
additions  are  made  to  the  same  sample  solution.  IDC-GSAM  was  found  to  be 


inappropriate  when  making  time  additions.  With  only  1%  random  noise, 
errors  as  high  as  30%  in  nQ  were  obtained.  The  inadequacy  of  the  IDC-GSAM 
to  differentiate  the  systematic  drift  from  the  1%  random  noise  is  an  arti¬ 
fact  of  the  experimental  design,  particularly,  the  perturbations  performed 
on  R.  Because  the  IDC-GSAM  volume-corrected  responses  can  be  successively 
increasing  thoughout  the  standard  additions,  the  noise  introduced  is 
subsequently  also  accumulative,  thereby  suggesting  that  the  IDC-GSAM  is 
very  sensitive  to  measurement  noise.  For  IDC-PGSAM,  the  sample  is  parti¬ 
tioned  yielding  additive  volume-corrected  responses  accumulating  for 


only  one  partition  and  reponses  are  reduced  back  to  the  initial  values 
before  a  series  of  standard  additions  are  made  to  a  new  partition.  Thus 
the  noise  introduced  is  no  longer  accumulative  for  the  period  of  the  analysis 
but  only  increasing  during  additions  to  a  partition.  Accordingly,  when 
time  additions  are  made,  the  1DC-PGSAM  is  the  method  of  choice. 

Three  studies  were  done  to  test  the  applicability  of  the  GSAM  drift 
correction  to  ICP-AES.  These  were:  no  drift  present  with  time  additions; 

10%  drift  in  the  A1  sensor  without  time  additions;  and  10%  drift  in  the 
A1  sensor  with  time  additions.  Table  VII  contains  the  final  results  for 
the  above  studies.  Spectral  interferences  occurring  were  the  well-documented 
interferences  o£  direct  spectral  overlap  of  As  (228.812  nm)  on  Cd  (228.802  nm) 
and  the  continuum  emission  of  A1  effecting  As  (193.696  nm)  (8).  A  previous 
study  demonstrated  the  applicability  to  GSAM  to  characterize  and  correct 
for  those  -.n. ‘terences  (3). 

The  results  of  this  study  indicate  that  when  time  additions  are  made, 

with  drift  present  or  not  present,  accurate  n  's  are  estimated.  With  the 

o 

occurrence  of  drift,  and  time  additions  are  not  included,  the  n  ’  s  are 

o 

largely  distorted.  Thus,  the  GSAM  can  be  used  to  accurately  correct  for 
instrumental  drift  in  ICP-AES  In  the  presence  of  sample  dependent  inter¬ 
ferences. 

Potential  Interferences  Study 

As  stated  in  the  Experimental  Section,  the  same  simulated  experiment 
used  in  the  time  study  with  no  drift  and  no  time  additions  (Table  I)  was 
used  with  one  minor  change  here.  That  is,  to  exemplify  the  simplicity  of 
this  extra  determination  of  potential  Interferences,  only  one  addition 
of  the  pseudo-analyte  was  made.  More  standard  addition  of  the  pseudo¬ 
analyte  may  be  performed  but  only  one  is  necessary  to  determine  the 
presences  of  potential  Interfering  analytes.  Even  though  in  the  Time 
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Study  sensor  3  represented  one  analyte,  this  sensor  can  also  be  the 
pseudo-sensor  obtained  by  averaging  the  responses  of  several  sensors. 

The  number  of  sensors  or  responses  constituting  the  pseudo-sensor  or 
pseudo-response  is  immaterial,  for  this  study,  as  long  as  the  average 
pseudo-response  equals  the  response  for  sensor  3  used  in  the  Time 
Study.  The  IDC-GSAM  could  have  been  used  here,  since  no  time  additions 
are  made;  no  statistical  differences  in  the  results  are  expected. 

Assuming  that  sensor  3  in  Table  VTII  is  the  average  of  many  sensors, 
then  the  GSAM  appears  to  correct  for  this  potential  interference.  k_ 

(last  row  of  k)  would  tell  the  analyst  that  sensor  two  could  have  a  sig¬ 
nificant  interference  pertaining  to  some  chemical  species  in  the  sample  if 
these  chemical  species  are  present.  k_  .  reveals  the  total  sensitivity 
of  all  sensors  making  up  the  pseudo-sensor.  The  nQ  associated  with  the 
pseudo-analyte  discloses  to  the  analyst  that  other  species  are  present  in 
the  sample  since  the  pseudo-nQis  significantly  different  than  zero.  Thus, 
further  qualitative  studies  of  the  sample  are  necessary  to  enable  the 

usage  of  the  appropriate  standard  addition  for  proper  quantitative  results 
of  the  secondary  analytes  composing  the  pseudo-analyte.  The  case  with  no 

extra  standard  addition  of  the  pseudo-analyte  was  also  investigated  with 

the  results  presented  in  Table  IX.  These  results  show  the  error  acquired 

in  the  determination  of  nQ  for  analyte  two.  Clearly,  when  possible,  it 

seems  best  to  make  one  extra  addition  of  all  conceivable  interfering  species 

and  average  the  respective  reponses  forming  the  pseudo-response  to  permit 

proper  correction  and  detection  of  the  potential  analytes. 
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CONCLUSION 

In  this  paper  an  extension  of  the  GSAM  was  presented  to  correct 
for  drift.  This  method  may  be  extended  to  higher  order  functions  of  drift 
provided  additions  of  equal  orders  are  made.  The  experimental  results  show 
that  the  IDC-PGSAM  may  be  applied  when  drift  is  uncertain.  With  the 
GSAM,  the  analyst  can  detect  and  correct  for  potential  interferents 
provided  at  least  one  standard  addition  containing  all  of  these 
interferents  is  made.  The  single  extra  addition  of  the  pseudo-analyte 
represents  a  useful  check  on  the  analysis.  If  the  extra  sensors  required 
for  the  check  on  the  analysis  are  available,  they  should  be  used  to  gain 

maximum  information  about  the  sample. 

The  incorporation  of  automatic  drift  detection  and  correction,  more 

sensors  than  analytes,  and  the  compensation  of  potential  interferents  in 
the  matrix  into  the  GSAM  represents  a  step  toward  the  development  of 
Intelligent  analytical  Instrumentation.  The  combination  of  computer  con¬ 
trolled  analytical  instrumentation  and  the  advanced  GSAM  as  the  chosen 
analytical  experimental  design  permits  instruments  to  perform 
self-correcting  multicomponent  analysis.  Studies  of  this  nature  are 
currently  underway  in  our  laboratory. 
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Table  I.  Experimental  Design  for  Simulated  Time  Study  (1DC-PGSAM) 


At 

At2 

- -p. 

partition  1 

1.0 

1.0 

1.0 

3.0 

1.0 

5.0 

partition  2 

5.0 

25.0 

1.0 

11.0 

1.0 

13.0 

partition  3 

9.0 

81.0 

1.0 

19.0 

1.0 

21.0 

K(without  drift) 


K(with  drift) 


12  3  1 


1 

1.0 

0.5 

0.0 

1.0 

2 

0.0 

1.0 

0.5 

0.0 

3 

0.0 

0.25 

1.0 

0.0 

t 

0.0 

0.0 

0.0 

0.02 

t2 

0.0 

0.0 

0.0 

0.0 

2 

3 

1 

2 

3 

0.5 

0.0 

1.0 

-0.5114 

0.2857 

1.0 

0.5 

0.0 

1.143 

-0.5714 

0.25 

1.0 

0.0 

-0.2857 

1.143 

0.0  0.0 

0.0  0.02 


1.0  0.5  1.0 


Table  II.  Experimental  Design  for  Al-As-Cd  Time  Study 
(IDC-PGSAM) 


AN (Mg) 


Al 

As 

Cd 

additions 

801.3 

0.0 

0.0 

of 

801.3 

0.0 

0.0 

Al 

801.3 

0.0 

0.0 

additions 

0.0 

599.4 

0.0 

of 

0.0 

599.4 

0.0 

As 

0.0 

599.4 

0.0 

additions 

0.0 

0.0 

101.8 

of 

0.0 

0.0 

101.8 

td 

0.0 

0.0 

101.8 

At (ml  minute) 

nQ(yg) 

36 

Al 

2503. 

48 

As 

999.7 

60 

Cd 

20.14 

168 

70 

82 


300 

92 

104 


Table  III.  Final  Results  with  no  Time  Additions,  no  Drift, 
and  1%  Random  Noise 


1.001  (0.016) 
-0.002406(0.0055) 
-0.004081(0.018) 


0.5029(0.015) 
1.007  (0.022) 
0.2541(0.048) 


-0.002333(0.0087) 
0.5043  (0.019) 

1.016  (0.040) 


K-1(s) 


0.9986  (0.014) 
0.0006847(0.011) 
0.003782  (0.012) 


-0.5715(0.018) 
1.138  (0.053) 
-0.2894(0.066) 


0.2862(0.015) 
-0.5655(0.044) 
1.130  (0.062) 


calcd.  n  (s) 
o 


relative  error  %c 


1  1.003  (0.024) 

2  0.4909(0.053) 

3  0.9903(0.058) 


0.300 

1.82 

0.970 


relative  error  =  100  x 


|  true  -  , al.  d. ' 


1 


Table  IV.  Final  Results  with  Time  Additions,  No  Drift, 
and  1%  Random  Noise 


1 

1.006  (0.012) 
0.001180  (0.0062) 
0.001514  (0.019) 
-0.007131  (0.012) 
0.0006920(0.0013) 


0.9946  (0.011) 
-0.0005412(0.010) 
-0.001299  (0.021) 


K(s) 

2 

0.4993  (0.016) 

1.003  (0.01V) 

0.2499  (0.049) 

0.01114  (0.021) 
0.001603(0.0027) 


K  1 2 3(s) 


-0.5662(0.016) 
1.142  (0.040) 
-0.2858(0.060) 


3 

0.001755  (0.0085) 
0.5043  (0.012) 

1.003  (0.046) 

0.006275  (0.0081) 
0.0005720(0.0016) 


0.2832(0.012) 
-0.5753(0.030) 
1.144  (0.066) 


calcd.  n  (s)  relative  error  % 
o  _ 


1  0.9889(0.023) 

2  0.5041(0.042) 

3  0.9955(0.068) 


1.11 

0.820 

0.450 


Table  V  Final  Results  with  No  Time  Additions,  2% 
Quadratic  Drift  in  Sensor  3,  2%  Linear 
Drift  in  Sensor  1,  and  1%  Random  Noise 


K(s) 


1 


3 


1  1.023  (0. 

2  0.04745(0. 

3  0.07381(0. 


017) 

0. 

0045) 

0. 

0054) 

0. 

4967(0.0080) 

9976(0.014) 

2461(0.0055) 


0.5799(0.0058) 
0.8229(0.010) 
1.798  (0.025) 


K~1 2 3(s) 


1  0.9875  (0.015)  -0.5456(0.015) 

2  -0.01520(0.0046)  1.139  (0.019) 

3  -0.03845(0.0033)  -0.1335(0.0055) 


0.2179(0.011) 

-0.5207(0.015) 

0.6185(0.0099) 


calcd.  n  (s)  relative  error  % 
o  _ 


1  0.9198(0.022) 

2  0.7212(0.021) 

3  0.3365(0.018) 


8.02 

44.2 

66.3 


Table  VI 


Final  Results  with  Time  Additions,  2%  Quadratic  Drift 
in  Sensor  3,  2%  Linear  Drift  in  Sensor  1,  and  1% 
Random  Noise 


K(s) 

12  3 


1  1.005  (0.014) 

2  0.002146  (0.0060) 

3  0.0001118(0.027) 

t  0.01990  (0.0032) 

t2  0.0003465(0.0011) 


0.5001(0.013) 

0.9990(0.016) 

0.2593(0.037) 

-0.004744(0.0069) 

0.001060(0.0054) 


0.0004186(0.0099) 
0.5009  (0.013) 
1.018  (0.061) 
0.008284  (0.0092) 
0.01905  (0.0022) 


^(s) 


1  0.9964  (0.012)  -0.5723(0.017) 

2  -0.002449  (0.017)  1.151  (0.041) 

3  0.0003758(0.030)  -0.2962(0.057) 


0.2817(0.016) 
-0.5677(0.034) 
1.132  (0.084) 


calcd.  n  (s) 
o 


relative  error  Z 


1  0.9945(0.022) 

2  0.4917(0.049) 

3  0.9893(0.082) 


0.550 

1.66 

1.07 


Table  VII  Final  Results  for  Al-As-Cd  Time  Study 


no  drift  2429. 

with  time 
correction 

10%  drift  2125. 

without  time 
correction 

10%  drift  2429. 

with  time 
correction 


calcd.  nQUg 


1020. 


1028. 


1020. 


20.18 


19.63 


20.18 


relative  error  % 


A1  As  Cd 


2.95  2.03  0.199 


15.1  2.83  2.53 


2.96  2.03  0.199 


Table  VIII.  Final  Results  of  Unknown  Interference  Study 


K(s~> _ 

1  0.9987  (0.017)  0.5010(0.0070)  -0.0005370(0.0053) 

2  0.002160(0.020)  1.007  (0.025)  0.4983  (0.0092) 

3  0.005021(0.055)  0.2415(0.074)  1.007  (0.043) 

calcd.  n  (s) 
o 


1  0.9937(0.063) 

2  0.4990(0.096) 

3  0.9958(0.063) 


relative  error  % 

0.630 

0.200 


Table  IX.  Final  Results  of  Unknown  Interference 
Study  without  Pseudo-sensor 


K(s) 


1  2 

0.9981  (0.015)  0.4974(0.011) 

-0.0002085(0.0052)  1.006  (0.015) 


calcd.  n  (s) 
o 


relative  error  % 


1.002  (0.022)  0.200 

0.7457(0.023)  4.91 
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